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expansion formalism. The essential attention has been given to such relativistic 
problem as a deuteron wave function in a moving frame and transformation of spin 
states due to Wigner rotation. Parameterization of the nucleon-nucleon t-matrix 



^ I has been used to take the off-energy shell effects into account. The vector, Ay, and 



tensor, Ayy, analyzing powers of the deuteron have been calculated at two deuteron 
kinetic energies: 395 MeV and 1200 MeV. The obtained results are compared with 



> 

T^ij- I the experimental data. 
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1 Introduction 

The study of the deuteron-proton elastic scattering has a longtime story both from the 
experimental [T]- [B] and theoretical [7|- [S] (and refs. therein) points of view. This process 
is the simplest example of the hadron nucleus collision that is why the interest to this 
reaction is justified. A number of experiments on deuteron- nucleon elastic scattering is 
aimed at getting some information about the deuteron wave function as well as nucleon- 
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nucleon amplitudes from Nd scattering observables. Moreover, the study of the reaction 
mechanisms, investigations of the few-body scattering dynamics are also very important. 

A good theoretical description of the deuteron-nucleon process is obtained for low 
energies, where the multiple scattering formalism based on the solution of the Faddeev 
equations has been applied to this problem [7j. Nevertheless, this method requires the 
realistic nucleon- nucleon potential but the well known A^A^ potentials such as Bonn [TU], 
CD-Bonn [IT], Paris [12], Nijmigen [13], AV18 [H], can be used only up to about 300- 
350 MeV A^A^ laboratory kinetic energy. This fact indicates a necessity to develop new 
relativistic nucleon-nucleon potentials which describe existing A^A^ data up to 1 GeV and 
higher. Many efforts to get these kind of A^A^ potentials have been undertaken |15j-[T9] 
but these potentials describe very well only low A^A^ phase shifts. It is possible in the 
future more acceptable relativistic A^A^ potential will be developed what allows to extend 
energy region where Faddeev equations formalism can be applied. 

At present there are experimental data on the dp elastic scattering in a wide range, 
where the deuteron kinetic energy is above 400 MeV [1]- [6] (and refs. therein). Already 
from a few hundred MeV of nucleon-nucleon interaction energy the relativistic effects 
begin to play a very important role. A relativistic quantum theory including relativistic 
kinematics, Lorentz transformation and the relativistic spin theory should be applied to 
consider a few-body scattering. The basic principles of this theory were formulated in Ref. 
[20] where results of the works of Dirac [21], Wigner [22j, Foldy [23], etc., were developed. 
In Ref. [21] review devoted to the relativistic dynamics is presented. Recently relativistic 
kinematics, boost effects and Wigner rotation have been included into the three-nucleon 
Faddeev equations [25], [26] . 

The present paper has realized the semi-phenomenological method based on the pa- 
rameterizations of the deuteron wave function (DWF) and nucleon-nucleon vertices. At 
intermediate energies the relative momentum of the two nucleons in the deuteron is rather 
large and, in general, the use of a relativistic deuteron wave function is necessary. In Refs. 
|27] . [28] these DWF's were obtained by solution of the Blankenbecler-Sugar equation. 
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But in this approach one of the nucleons is on-energy shell that limits the employment 
of these DWF's. The two nucleons in the deuteron are considered as the off-energy shell 
particles in Ref. [29], where the deuteron wave function is the result of the Bethe-Salpeter 
equation solution. The importance of the relativistic DWF application to describe the 
deuteron-proton backward elastic scattering was demonstrated in this paper. 

However, at the energies considered in the present paper it is possible to use non- 
relativistic DWF such as Paris [30], Bonn [10], CD-Bonn [11] due to transition into the 
deuteron Breit frame. It allows one to decrease the two nucleon relative momentum both 
in the initial and final deuterons. In this paper the CD-Bonn DWF has been used as 
a more modern one. It should be noted that the Lorentz transformation and relativis- 
tic spin theory have been applied to describe the transformation from the deuteron rest 
frame to the Breit system. As a result, the deuteron wave function has a more compli- 
cated spin structure, than usual nonrelativistic DWF and depends on the two variables. 
The analogous result was obtained in Refs. [3T]-[33] for the DWF on the light front. 

The model offered in Ref. [8] to investigate the cip-elastic scattering at intermediate 
energies, has been developed further in the present paper. First of all, the deuteron wave 
function transformation has been done. Secondly, the phase shift analysis data are not 
used to describe the nucleon-nucleon interaction in comparison with Ref. [8], where the 
energy dependence in A^A^- vertices important for correct integrations, is lost. Besides, the 
phase shift analysis gives information about on-energy shell nucleons. But in our method 
the off-energy shell effects have to be taken into consideration and it can be done by 
using some parameterization describing the nucleon-nucleon interaction. Thus, the A^A^ 
t -matrix constructed by Love and Franey [31] has been used in our calculations. 

The presented method has been applied to calculate two simplest polarization observ- 
ables: vector. Ay, and tensor, Ayy, analyzing powers at two deuteron kinetic energies: 395 
MeV and 1200 MeV. The results of these calculations are compared with the experimental 
data obtained in Saclay [1], [2] and Argonne [3|. It should be noted that these observ- 
ables and energies are examples of the presented method employment. In the future this 
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method is planned to describe other polarization observables at deuteron kinetic energies 
between 400 and 1200 MeV. 

The paper is organized as follows. Section 2 gives the general formalism: the transfor- 
mation of the deuteron wave function from the rest frame to a moving deuteron system, 
is considered; the description of the nucleon-nucleon interaction is presented; the calcula- 
tions of the scattering amplitude terms are performed. Section 3 defines the polarization 
observables, Ay and Ayy. The results of the calculations are discussed in Sect. 4. The 
conclusions are contained in Sect. 5. 

2 General formalism 

According to the three-body collision theory, the amplitude of the deuteron-proton elastic 
scattering J' is defined by the matrix element of the transition operator Uu: 

Udp^dp = 5iEd + Ep-E'a-E'p)J=<l{23)\[l-Pi2-Pi3]Uu\li23)> . (1) 

Here the state |1(23) > corresponds to the configuration, when nucleons 2 and 3 form 
the deuteron state and nucleon 1 is free. Appearance of the permutation operators for 
two nucleons Pij reflects the fact that the initial and final states are antisymmetric due 
to exchange of the two particles. 

The transition operators for rearrangement scattering are defined by the Alt-Grassberger- 
Sandhas equations [35], [36] 

Uu = t2goU2i + hgoU-ii 

U21 = 9o^ + hgoUn + t3goU3i (2) 
U31 = + hgoUn + t2goU2i , 

where ti = t(2,3),etc., is the t-matrix of the two-nucleon interaction and qq is the free 
three-particle propagator. The indices ij for the transition operators Uij denote free 
particles i and j in the final and initial states, respectively. 

Iterating these equations up to the ti second order terms we can present the reaction 
amplitude as a sum of the following three contributions: one nucleon exchange, single 
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scattering and double scattering terms, - 

Jdp-^dp — JoNE + JsS + JdS 

JoNE = -2<l(23)|Pi2(?o'i|l(23)> 
Jss = 2<l(23)|ini(23)> (3) 
Jns = 2<l(23)|tf"^^otr™|l(23)> , 

where we have introduced notations for antisymmetrized operators ^2^™ = [1 ~ Pi?]t2 and 
fy'^ = [1 — Pi2]H- Before we consider each of these elements in detail, we must determine 
the reference frame. 

Breit system 

The most popular deuteron wave functions are well known only in the nonrelativistic 
region. In this connection it is very important to choose the frame, where the application 
of the nonrelativistic deuteron wave function is possible. The frame which minimizes the 
relative momenta of the nucleons in the both deuteron wave functions, is the deuteron 
Breit system. In this frame the deuterons move in opposite directions with equal momenta 
(fig.l). As a consequence, the energies of the initial and final deuterons (and protons) are 
equal to each other. 

{pQ) = -Q' ■ (4) 
We take the orthonormal basis 

z^j^ — z^ = -Q: X = —-— = p + Q, y = zxx=pxQ , (5) 
\p — p\ \P + p\ 

where ^-axis goes along the transfer momentum and x-axis — along the average momen- 
tum of the initial and final protons. As usual, the y-axis is normal to the scattering 
plane. 
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Deuteron wave function 



In order to get the wave function of the moving deuteron, let us consider an arbitrary 



frame where two nucleons have momenta pi, p2 and energies Ei = Jrnj^ E- 



m'lf +P2^, respectively. Here mjsf is the nucleon mass. These momenta are related with 
those in the center-of-mass (cm.) by Lorenz transformation 

Liu)p, = {E*,p) 

L{u)p2 = iE*,-p) (6) 

with velocity 

P1+P2 
El + E2 

Here the cm. energy of one of the nucleons E* is related with Mandelstam variable s by 

^* = v^/2 . (8) 
Thus, we can introduce new variables Q and k which can be expressed through pi and p2 

Q =Pl+P2 

iE2 + E*)pi-iEi + E*)p2 
E1 + E2 + 2E* ■ 

Then a two- nucleon state in the (^1,^2) system is connected with a two- nucleon state in 

the center-of-mass by the relation 

\pi,P2>=J-'^\pi,P2)Wi/2{pi,u)Wy2{P2,u)\k,Q> , (10) 

where k is the relative momentum of the two nucleons in the cm. The normalization factor 
provides conditions of orthonormality and completeness and is defined by the Jacobian of 
the transformations 



The Wigner rotation operator in the spin space of the i-th nucleon has the standard form 

Wi/2{pi,u) = exp{-iui{niai)/2} = cos{uJi/2)[l - i{niai)tg{uJi/2)] (12) 
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with the rotation axis 

Ui = —r (13) 

\u X pi\ 

Ref.pO] shows that under a pure Lorenz transformation 

L{u')P = Q (14) 

the wave function of the bound state transforms in the same way as the state of a single 
particle 



W{L^,)\P>-- 



. ^ W,{Q,u')\Q> , (15) 



where the Wi is the Wigner rotation operator for spin 1 particle 

Wi{Q,u') = exp[-iuj{n'S)} (16) 

with the rotation axis 

^, u' xQ 

n = — . 17 

\u' xQ\ 

The deuteron wave function in the rest frame of the deuteron depends on only one 

variable k, which is the relative momentum of the outgoing proton and neutron 

1 f w(k) ^ ^ ] 

< mpmr,\nd\Md >= -j= < mprrinl { u{k) + -^[3{aik){a2k) - (^1^2)] \ \Md > ,(18) 



/47r 

where u{k) and w{k) describe the S and D components of the deuteron wave function [10], 
[TT] . [30] . k is the unit vector in k direction and A4d,fnp,fnn are spinors of the deuteron, 
proton, and neutron, respectively. 

Obviously, the deuteron rest frame corresponds to the center-of-mass of the two out- 
going nucleons, P = Pn + Pp = 0. Then the proton -neutron relative momentum k in 
Eq. lfTSi) is related with the nucleon momenta pi and p2 in the arbitrary frame by Eq.(l9]). 
Using the formulae (p!)- (fT8|) we get the expression for the deuteron wave function in the 
frame, where the deuteron is moving and its momentum is Q 



< Pi P2,mim2\fld\Q,Md > 



Ma E*{E^ + E2) uji 002 
< A;Q,m;m2|[l -2(frai)t^y][l + i(nCT2)t^Y]fid|0,A^d > . 



The Wigner rotation is performed around the normal to the (kQ) plane 

Pixp2 kxQ 
n = — — = — — (20) 

\Pixp2\ \kxQ\ 

on the rotation angles, which are defined by the relations |37j : 

tg^ = , (21) 

^ <^2 \P\ X P2I 



2 mjv(^i + ^2 + 2E*) + 2E*(E2 + E*) 
The expression (1191) can be presented through six terms: 

<PxP2-,mxm2\^d\Q-,Md> = gi{k,Q) + g2ik,Q){ain){a2n) + g3{k,Q){aia2) + 

+ 9dk,Q){aik){a2k)+g5{k,Q)[{ai + a2)n]+ (22) 
+ g6{k,Q)[iaik){a2n x k) + (aiu x k){a2k)]. 



Solving a relativistic equation one can obtain functions gi as it was done in Ref. 
on the light front. But in this paper a usual nonrelativistic DWF is taken as input. 
Therefore, functions gi are defined as linear combinations of u and w {S- and D-waves). 

This way we have got the relation between the deuteron wave function in two systems: 
the deuteron rest frame and the arbitrary moving system. Note, the Wigner rotation 
corrections are small at the energies under consideration, but momenta transformations 
are very important since they give us an opportunity to use the nonrelativistic deuteron 
wave functions in a rather wide energy region. 



One nucleon exchange scattering 

The first contribution into the rfp-elastic scattering amplitude in Eq.([3]) is the one 
nucleon exchange (ONE) term. The corresponding diagram is presented in Fig. la. Ap- 
plying the definitions of the wave function of a moving deuteron and three- nucleon free 
propagator, we can write ONE amplitude in the following form: 

JoNE = -2 1(23) < pm'r'; -QM'aO\n[{23)Pu (23) 
{Ed + Ep-ki- k2 - + ie)VLd{2?>)\QMd^;pmT >^2z) , 
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where m, m' are the spin projections of the initial and final protons, r, r' are their 
isospin projections, respectively. The kinetic-energy operator has the standard definition. 



Ki\pi >= ym%+pi^\pi > . 

After a straightforward calculation we have the expression for the ONE amplitude 



JoNE = -^{Ed - Ep ~ + p^ - Q^) ■ (24) 

< 0m'; -QM'a\nli{23)[l + {aia2)]nd{23)\QMd;pm > , 

where the definition of the permutation operator in spin space Pi2(cr) = |[1 + {cTicr2)] 
has been applied. The two nucleons relative momentum in the deuteron rest frame is the 
argument of the deuteron wave function in Eq.l fTSj) . Using Eq.Q we can express these 
variables for the initial Pq and final ff^ deuterons, correspondingly, 

En + E* 



Po 



p + Q 
p + Q 



1 + 



1 - 



Ep + E^ + E* 
En + E* 

Ep + En + E* 



(25) 



Here En = y mj^ + p"^ — Q'^ and E* = y {Ep + EnY — are the struck neutron energy 
in the moving deuteron frame and rest deuteron frame, respectively. Note, that \pq\ = \j7q\. 
The Wigner rotations are performed around y-axis on the equal angles both for the initial 
and final deuterons. Using Eq. (l2Tl) we get 



0J\ 



\P^Q\ 



tg 



UJ2 



tg- 



2 mN{Ep + En + 2E*) + 2E*{Ep + E*) 
IpxQI 



(26) 



2 ""2 mN{,Ep + En + 2E*) + 2E*{En + E*) ' 
where the angles wi, uj[ are related with the projectile and outgoing protons, respectively, 
while U2, uj'2 correspond to the struck neutron (Fig. la). 



Single scattering 

The other term in the dp-elastic scattering amplitude Eq.(l3]) is the single scattering 
one. The corresponding diagram is presented in Fig. (lb). Following the standard proce- 
dure and applying the definition of the moving deuteron wave function Eq.( |T9l) . we get 
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the expression for the single scattering amphtude: 



Jss = J dq' <-QM'MW'm",-Q-q'm',> (27) 
< p'm, -Q - q'\-tl2 + 2^i2\wn^ Q - q'm^ X q'm", Q - q'm2\ild\QMd > ■ 

The relative momenta of two nucleons for the initial and final deuterons are 

Po = q +Q 



E2 + E'3 + 2E'* 



where the nucleons energies E2, E3, E'^ in the reference frame are defined by the standard 
manner (Fig. lb) 



E2 = ^m% + q'^, Es = ^m% + iQ~q'y , E', = ^ m% + (Q + q (29) 

and these energies in the center-of-mass of the two nucleons forming the initial and final 
deuterons are equal, correspondingly, to 



E* = -^{E2 + E,f - g2 , E'* = -^{E2 + E',Y - g2 . (30) 

The Wigner rotation axis goes along the normal to {q'Q) plane 

^ q' xQ 

n = — 31 

Ig'xQI 



and Wigner rotation angles are defined in accordance with Eg. (1211) for the initial 

, ^ _ \q' X Q\ (^2) 

^2 mN{E2 + E3 + 2E*) + 2E*{E2 + E*) ^ ' 



2 mNiE2 + E3 + 2E*)+2E*{E3 + E*) 

and final deuterons 

, ^ W X Q\ (00) 

^2 mN{E2 + E^ + 2E'*) + 2E'*{E2 + E'*) ^ ' 

, ^ _ 1^' X Ql 



2 mN{E2 + E'3 + 2E'*) + 2E'*{E^ + E'- 
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Further we have to give essential attention to the two nucleon t*-matrix in Eg. (1271) . 
Here index i denotes the isospin state of the nucleon-nucleon system. This t-matrix is 
the half-off-energy shell one. Since we assume to use the formalism presented here to 
investigate the processes at rather high energies, we can not define the high-energy A^A^ 
t-matrix solving the Lippmann-Schwinger equation with any known potential. In this 
situation it is possible to apply some parameterization as in Refs. |3S], [3S] to study the 
dp breakup in the GeV region. 



Nucleon-nucleon t-matrix 

In order to describe the nucleon- nucleon interaction in a wide energy region, we use 
the Love and Franey parameterization [31]. In this model the on-shell NN t-matrix is 
defined in the center-of-mass 

< K'm[m'2\t\K,mim2 > = < K'm[m2\A + B{aiN*){a2N*) + C{ai + ■ N* + 

D{aiq*){a2q*) + F{aiQ*){a2Q*)\Kni^ni2 > , 

where the orthonormal basis is combinations of the nucleons relative momenta in the 
initial k and final k' states 

q* = T^^ , Q* = j^-^ , N'' = —— . (34) 

\K — K \ \K + K'\ \K X k'\ 

We can not directly use this definition of the t-matrix in our calculation. First of all, 
because this t-matrix is defined in the cm., while we need that in the deuteron Breit frame. 
Second, we need the off-energy shell t-matrix, but the Love and Franey parameterization 
gives us the on-energy shell matrix. However, the both problems can be solved. 

Using formulae f|T0|) - f|T3|) we can connect expressions for the NN t-matrix in the cm. 
with the same matrix in the reference frame by the following relation: 

< pp;^;m'm'^\t\pp3;mm3 >= J\f < K'm^m2|Vr//2(p') ^1/2(^3) ^ Wi/2{p)Wi/2{p3)\n'mim2 > 

Here momentum notations correspond to the ones given in Fig. lb. The Wigner rotations 
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of the initial and final states are performed around jl and jl' -axes, respectively: 

A* = — ^5 A* = "FTT ^ • (36) 

We introduce the following notations for the Wigner angles: xi^ x'l correspond to the 
projectile and outgoing protons, while Xs? Xs ^ire related with the nucleons from the 
initial and final deuterons, correspondingly. These angles can be expressed through their 
tangents using Eq. fl2T|) : 



. Xi |P X PsI 



2 rriN (Ep + Es + 2El^) + 2^1*3 {E^ + ^1*3) 



^ X3 ^ \p X pz\ 



2 m;v(^p + ^3 + 2^1*3) + 2^1*3(^3 + ^1* 



13J 



, Xl P X P3 



2 (Ep + + 2^1*^) + 2^1*^ (Ep + £;r^) 



, Xs' \p' X ps'l 



2 m;v(i5p + i5^ + 2i5;i*^) + 2Ei*^(i5;^ + Ei*^) 

As it is shown in ref.[H] the spin rotation corrections to the A^A^ t- matrix are very 
small at the nucleon kinetic energy of 800 MeV . Therefore in our calculations we neglect 
the second term in Eq.( fT2il . Thus, only cosines of the Wigner angles are the Wigner 
operators in Eg. (1351) and, therefore, the A^A^ t- matrix in the reference frame has the 
same spin structure as that in the cm. 

The normalization factor in Eq. ( |35ll is related with the Jacobian of the transformations 
Eq. (fTTI) and equals to 



Eh{Ep + E3 



Et'iiEp + -E3) 



(38) 



^ 2EpE^ \ 2EpE'^ 

where E^^ and E'^^ are the initial and final cm. energies of the two nucleons taking part 
in the interaction, respectively: 

^13 = l^iEp + E^r-K^ , = yiEp + E^)^-K^ , (39) 

with a full momentum 

K=P+Q-^ . 
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In our calculations we need the half-off-energy shell and off-energy shell nucleon- 
nucleon t-matrices, while Eq. (!34|) gives us only the on-energy shell matrix. In order 
to reach this purpose, we use the result of the Ref. [40j which shows the off-energy shell 
t-matrix to be connected with the on-energy shell matrix by the relation 

t{Eonsheiu K; kk) = Ft{y/s; kk) , (40) 

where the kinematical factor F is defined as 

^ Eonshell ^s + 2E* ^S + 2E'* 



\fs Eonshell + Ep + E3 Eonshell + Ep + E'^ 

Here Eonsheii denotes the energy, which is determined by the following condition: 



(41) 



Eonshell = Ep + E^ = Ep + E'^ . (42) 



For the single scattering this condition will be satisfied at {Q(f) = 0. Then, E, 



onshell 



Ep + ,Jm% + Q^ + q'^. 

The argument of the on-energy shell t-matrix in Eq. (HO!) is the square root of the 
Mandelstam variable s: 

S = ^onsheZ« - • (43) 

Thus, using the scenario presented here we can extend the parameterized on-energy 
shell A^A^ t-matrix to the off- energy shell case. 

Double scattering 

The diagram of the double scattering is shown in Fig.(lc). The corresponding ampli- 
tude J'ds can be written as 

Jds = NN' I dq' I df <-QM'a\nl\q'm'^,Q-q' m'^> (44) 



')] [tUV^, k, k') + 4(v^, K k')] 

1 ^ 

\mm2mz>< -Q - q'' ■m2.,q'' ■mz\Vtd\Md> ■ 

13 



Ed + Ep-E[-E'2-E'^ + le 



As before, we define tlie relative momenta of the nucleons which form the initial and final 



deuterons, respectively, in accordance with Eq.([9]) 

where the center-of-mass energies of these nucleons are related with those in the reference 
frame by 



E* = -^{E2 + E,y-Q', E'* = {E', + E',y - . (46) 

Here the energies E2, E3 correspond to the nucleons in the initial deuteron. The expres- 
sions for them are identical to the ones in the single scattering case and defined by Eq. ( l29ll . 
The primed energies E[,E2,E'^ are related with the corresponding momenta (Fig.lc) in 
the standard manner 



E[ = ^m"^ + {p + Q -q' -q")\ E'^ = m"^ + [Q + q , E'^ = ^ + q .{47) 

Applying Eqs. fl20|) - fl2T]) to describe the transformation of the deuteron wave function 
due to the Wigner rotation, we can find the directions of the rotation axes for the initial 
and final deuterons, respectively: 

^ q' xQ ^, q" X Q 

n = — , n = — . 48 

Ig'xQI |g"xQ| 

Note, these axes do not coincide here as it was in the one nucleon exchange and single 

scattering cases. 

The rotation angles are expressed through their tangents, which are defined by Eq. (l32|) 
for the initial deuteron, and by the following relations — for the final deuteron: 

tq^ = l^"^^^^! (49) 

^2 m{E'2 + E'^ + 2E'*) + 2E'*{E!2 + E'*) ^ ' 

^ _ \q"^Q\ 
9 o 



2 m{E'2 + E'3 + 2E'*) +2E'*{E^ + E'*) 

In Eq. (jH]) we have two A^A^ t-matrices which correspond to two rescattering vertices. 
There are the off-energy shell matrices defined in the previous subsection. We get a/si3 
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and \J s'i2 as the arguments of the first and second t-matrices, correspondingly. The cm. 
relative momenta of the colliding and scattering nucleons in the first vertex are related 
with the nucleon momenta in the deuteron Breit frame by the following expressions: 

(p + Q){E^ + El^) + q'{E, + El,) 



k 
k' 



Q- 



Ej) + E3 + 2Ef- 



(50) 



q"-{p + Q-q" ' 



13 

rpl I I?'* 



The arguments of the second NN t-matrix are 



q" - {p + Q-q' 
-Q 



E[ + E', + 2E 



E2 + El 



E[ + E2 + '2EI2 
{p + Q){E'2 + E'*2) + q"{E, + E[*2) 



(51) 



Er, 



E'2 + 2E'{2 



where the two-nucleons center-of-mass energies can be expressed by means of the corre- 
sponding Mandelstam variables s: 



7-1* 
E12 



7-1/* 

-^12 



= ]^^{E, + E2y-[p + Q-q'r 



s'i2/2 = -V(^P 



M)2 



{p + Q-q") 



(52) 



E- 
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7-1/* 

-£^13 



S13/2 

^/2 



[E. + E^y-ip + Q-q') 
{E, + E',Y-{p + Q-q') 



The normalization factors and A^' are connected with the transformations of the 
nucleon-nucleon t-matrices as follows: 



-^1*3 (-^1 + E',, 



EhjEp + -^3 

^ 2EpE3 \ 2EiE', 



'Sl3 X X3 Xi X3 
—cos — cos COS COS — 

Ep + E3 2 2 2 2 



N' 



lEl2{Ei + E2) 



EUEp + E'2) V*12 (^1 ^2 ^ ^'2 

A c^Tp -n^i TT ^COS — COS — COS — COS — 

\ 2EpE2 Ep + E2 2 2 2 2 



(53) 



(54) 



y 2E1E2 

Each of these factors includes the normalization coefficient M related with the Jacobians 
of the transformations from the cm. to the reference frame, Eq. ffTTj) . the kinematical 
factor F defined in Eg. (1411) . and cosines of the Wigner rotation angles which follow from 
the Wigner rotation operators, Eg. (1121) . The expressions for the tangents of the Wigner 
angles are given in Appendix. 
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The three nucleon free propagator in Eq.(jH]) can be decomposed on the two terms 
using the well known formula: 

, , . = V \ , ,- 

Ed + Ep — E'l — E2 — E'^ + ie Ed + Ep — E'l — E'2 — E'^ 

(55) 

ni5{Ed + Ep-E[~E'2-E'^) . 

The principal value part is often neglected to simplify further calculations. However this 
method is efficient only to evaluate the double scattering contribution. In order to get 
a more precise result, it is very important to take into account the both terms of the 
three nucleon propagator . In our calculations we perform all integrations without any 
simplifications by using standard CERNLIB programs. The difference between the results 
obtained with and without the principal value part of the free propagator, is discussed 
below. 

3 Polarization observables 

In the previous section the amplitude of the deuteron- proton elastic scattering has been 
considered in detail. We have shown how three contributions contained in this amplitude 
can be expressed in terms of the deuteron wave function and nucleon-nucleon t-matrices 
defined in the deuteron Breit frame. 

On the other hand, as in ref. [8j|, the dp elastic scattering amplitude Jdp^dp can be 
decomposed in spin-1 operators acting in the deuteron spin space and spin-1/2 operators 
cTj acting in the proton spin space: 

Jdp^dp = <rn'M'd\h + f2{Sy) + hQ^^, + fiQyy + h{ax){Sx) + 

f6{ax)Q^y + fri^y) + f8iay){Sy) + h{ay)Q^^, + (56) 

floi^y)Qyy + fll{^z){Sz) + fl2{^z)Qy:,\mMd > ■ 

Due to the time-reversal and parity invariance, we have only 12 linearly independent 
amplitudes fi. These amplitudes can be calculated by using the technique presented 
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above. The operators Uj are the Pauh matrices, while Si and Qij are the spin-1 vector 
and quadrupole operators, respectively: 

1 2 - 

Qij — —{SiSj + SjSi) — '^^ij^ 1 Qxx + Qyy + Qzz — . (57) 

This paper considers only two polarization observables: the vector analyzing power 
of the deuteron Ay and tensor analyzing power Ayy which are defined by the standard 
manner: 

_ TrjJSyJ^) TrjJQyyJ^) 

The squared dp-amplitude summarized over all the spin projections can be presented 
through amplitudes /j as 

rr(Jjt) ^ 6(A2 + /,^)+4(/| + /52 + /| + A\) + ^(/| + /| + /| + A^o)+ (59) 

(/e + /12) - \Mhfl + hKo) ■ 



Using the standard technique one can also obtain traces to define the vector 
Tr{JSyJ^) = Re 



8(/i/2* + Ms) - ^(/2/3 + /s/g) + (60) 



^(/2/: + /s/ro) + 2(/5/6 + fnfu) 



and tensor analyzing powers 



TriJQyyJ^) = lif.+fi) + lifs+f^)-lif! + a-lif! + fu)- 

lifi + fu) + iReifif: + fiKo) - \Re{hn + Ms) + (61) 
^i?e(/3/: + /9/ro) ■ 

In general, we have to apply the spin transformation technique to relate the laboratory 
observables to the Breit frame ones. However, we have considered only those observables 
which are defined by the polarization along the normal to the scattering plane. In this 
case the definitions of the polarization observables are identical both — for the laboratory 
and Breit frames. It is not the case for other observables, for example, A^x-, A^z, etc. 
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4 Results and discussions 



The results of our calculations are presented in Figs. (2-5). We have performed our in- 
vestigations at two deuteron kinetic energies: 395 MeV and 1200 MeV. This choice was 
motivated by two reasons. First, there are sets of experimental data at these energies 
[T]-[3]. Second, these energies have demonstrated the range, where the method presented 
here can be applied. 

The dashed curves in Figs. (2-5) correspond to the calculations which include only 
the one nucleon exchange and single scattering (ONE+SS) contributions. The solid line 
is the result which also takes the double scattering (DS) into account. All the calculations 
have been carried out with the CD Bonn deuteron wave function [11]. We have obtained 
a good agreement of the theoretical curves with the experimental data at the deuteron 
kinetic energy of 395 MeV. Note, the contribution of the double scattering is not very big 
at this energy, especially for Ayy. 

The other situation has been revealed for the higher energy. An agreement between 
the obtained curves and experimental data at the deuteron energy of 1200 MeV is not 
so obvious. The vector analyzing power is very well described up to the angle of 90*^ 
cm., while the behaviour of the tensor analyzing power is reproduced only qualitatively 
in a wide angular range. Here the single scattering contribution dominates only up to 
6* = 30^, whereas the both polarization observables are very sensitive at this energy to 
the double scattering at larger angles. Especially the behaviour of the Ayy changes very 
strongly due to the DS term having been taken into account. On the contrary, at the 
deuteron energy of 395 MeV the double scattering contribution is practically insignificant 
for the tensor analyzing power but for the vector analyzing power it is up to 50%. 

The results of the calculations with and without the principal value part of the three- 
nucleon free propagator are given in Figs. (6-9). The curves corresponding to the predic- 
tions for the deuteron energy of 395 MeV, are practically undistinguished. It is a rather 
clear result, since the double scattering contribution is not so large at this energy. 

At 1200 MeV the contribution of the principal value part of the free propagator is more 
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significant. Note, that the full calculation describes the vector analyzing power better than 
the simplified method, especially at the 40*^ — 80° cm. angles, where the double scattering 
is very important. But it is not the case for the tensor analyzing power. Here the full 
calculation gives the smaller absolute value as compared with the experimental data in the 
range, where Ayy has a maximum or minimum , while the result obtained by the reduced 
method is in a good agreement with the data. However, it is necessary to keep in mind, 
that this shortage is observed even in the range, where the single scattering dominates 
{6* < 40°). From the above it is possible to conclude, that the reason of this deficiency is 
not related with the three nucleon free propagator in the double scattering term. 

It should be noted that these observables are also sensitive to the nucleon-nucleon 
amplitudes. Unfortunately, to describe NN interactions we have used the model offered 
twenty years ago, when the nucleon-nucleon phase shift analysis was incomplete. Besides, 
up to now we have no reliable data of the np amplitudes at the nucleon kinetic energy 
higher than 1.3 GeV [H] and that is very important to study reactions with light nuclei at 
intermediate energies. Maybe, that the qualitative disagreement between the theoretical 
predictions and experimental data at the deuteron energy of 1200 MeV is caused by 
incomplete A^A^ interaction description. In this connection the modern parameterization 
of the nucleon-nucleon t-matrix in a wide energy range is needed. 

5 Conclusion 

In this paper we have presented a method to calculate the amplitude of the deuteron- 
proton elastic scattering at intermediate energies. Special attention was given to the 
questions connected with the relativistic effects. The transformation of the deuteron 
wave function in the rest frame to a moving system was performed, that allowed us to 
use the nonrelativistic DWF at rather high energies. In order to describe nucleon-nucleon 
interactions in a wide energy range, we have used parameterization of the A^A^ t-matrix. 
The spin transformation technique has been also applied to relate this t-matrix given in 
the cm. to that in the reference frame. 
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Using the method presented here we have managed to describe the experimental data 
on the vector, Ay, and tensor, Ayy, analyzing powers at two energies. We got a good 
agreement between our predictions and experimental data at the deuteron kinetic energy 
of 395 MeV and acceptable agreement — at 1200 MeV. We believe that the difference 
between the obtained predictions and data is related with the shortcomings of the nucleon- 
nucleon parameterization. It is possible that the deuteron-proton scattering description 
can be improved, if the modern NN model is used. 

In the future it is planned to apply the approach developed in this paper to describe 
experimental data at other energies. In particular, the data obtained at the Nuclotron 
(Dubna) at 880 MeV of the deuteron energy |42] can be described in this formalism. 

Acknowledgements: The author is grateful to Dr. V.P. Ladygin for fruitful discussions. 
This work has been supported by the Russian Foundation for Basic Research under grant 
A^^ 07-02-00 102a. 

Appendix 

Applying Eq. f|T2|) we can obtain the expressions for the Wigner angles. For the first 
nucleon-nucleon scattering we have 



, X 
to— 

^2 


|P X 


Psl 


m(Ep + ^3 + 2^1*3) 


+ "^El^iEp + -^1*3) 






PsI 




+ 2^1*3(^3 + ^1*3) 


tg^ 

^ 2 


\Pl X 


q"\ 


m(Ei + EI, + 2E'^,) 


+ 2EUEi + E[%) 


tg^ 

^ 2 


\Pl X 




m{E, + E:, + 2E[%) 


+ 2E[%{E', + E[%) 



where Xs ^ire the angles corresponding to the beam proton and struck nucleon from the 
deuteron, respectively, and Xi, x's are the angles corresponding to the scattered nucleons 
(Fig. Ic). 
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The analogous expressions are obtained for the second scattering vertex: 




m{E^ + E2 + 2Eh) + 2^1*2(^1 + ^1*2) 



m{Ei + E2 + 2EI2) + '2EUE2 + £^1*2) 

W X P2\ 



m{Ep + E!, + 2E[%) + 2K2(^p + K2) 
W X p'sl 



Here tpi, ip2 correspond to the coUiding nucleons, while 93, correspond to the scattered 
proton and nucleon from the deuteron, respectively. 
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Figure captions 



Fig.l The diagrams included into consideration, (a) The one nucleon exchange dia- 
gram, (b) The single scattering diagram, (c) The double scattering diagram. 

Fig. 2 The deuteron vector analyzing power at 395 MeV as a function of the cm. 
angle, the data are taken from ref. [1]. The dashed line corresponds to the ONE+SS, the 
solid line shows the ONE+SS+DS calculation. 

Fig. 3 The deuteron tensor analyzing power at 395 MeV as a function of the cm. 
angle, the data are taken from ref. [T]. The dashed line corresponds to the ONE+SS, the 
solid line shows the ONE+SS+DS calculation. 

Fig. 4 The deuteron vector analyzing power at 1200 MeV as a function of the cm. 
angle, the data are taken: (•) from ref. |2J, (A) from ref. [3j. The dashed line corresponds 
to the ONE+SS, the solid line shows the ONE+SS+DS calculation. 

Fig. 5 The deuteron tensor analyzing power at 1200 MeV as a function of the cm. 
angle, the data are taken: (•) from ref. [2j, (A) from ref. The dashed line corresponds 
to the ONE+SS, the solid line shows the ONE+SS+DS calculation. 

Fig. 6 The deuteron vector analyzing power at 395 MeV as a function of the cm. 
angle, the data are taken from ref. jlj. The solid and dashed lines correspond to the 
calculations with and without the principal value part of the propagator. 

Fig. 7 The deuteron tensor analyzing power at 395 MeV as a function of the cm. 
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angle, the data are taken from ref. [T]. The sohd and dashed hnes correspond to the 
calculations with and without the principal value part of the propagator. 

Fig. 8 The deuteron vector analyzing power at 1200 MeV as a function of the cm. 
angle, the data are taken: (•) from ref. (A) from ref. p]. The solid and dashed lines 
correspond to the calculations with and without the principal value part of the propagator. 

Fig. 9 The deuteron tensor analyzing power at 1200 MeV as a function of the cm. 
angle, the data are taken: (•) from ref. [2], (A) from ref. [3j. The solid and dashed lines 
correspond to the calculations with and without the principal value part of the propagator. 
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Figure 1: Fig.l 
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